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ACT I: Space Warps 




The Challenge of Interstellar Flight 








HBLtOPAUSe 


vOYAoem 



Voyager 1 mission: 

0.722 1 spacecraft launched 
in 1977 to study outer solar 
system and boundary with 
interstellar space. 

- After 30+ years, Voyager 1 
is currently at ~120 
Astronomical Units (AU) 

' - ^from the sun travelling at 
3:6 AU per year, 

I 

' IfN^yager 1 were on a 
-Hr^jectory h^ded to one of 
the Sun's ndcff^lf 
neighboring star systems. 
Alpha Centauri at 4.3 light 
years (or 271,931 AU), it 
would take ~75,000 years to 
traverse this distance at 3.6 
AU/year. 


Pictures courtesy of NASA, 

: i. nasa. 

http://voyager.jpl.nasa.gov/mission/interstellar.html 


DAEDELUS 






Project Daedelus 
sponsored by 
British 

Interplanetary 
Society in 1970's to 
develop robotic 
interstellar probe 
capable of reaching 
Barnard's star, at ~6 
light years away, in 
50 years. 

The resulting 
spacecraft was 
54,000t, 

92% fuel for fusion 
propulsion system. 

ISS is ~450t 



Picture courtesy of Adrian Mann, www.bisbos.com used with permission 
ISS picture courtesy of NASA (STS-135 fly-around) 
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Hyper-fast interstellar travel 


N'i^SA 


Is there a way within the framework of physics such that one 
could cross any given cosmic distance in an arbitrarily short 
period of time, while never locally breaking the speed of light 
(ll^h commandment)? 

SPACEWARPS 

(inflation) 


WORMHOLES 

(shortcuts) 


Images courtesy NASA 





Inflation: Alcubierre Metric^ 







Warp Metric: 



1 . 


Alcubierre, M., "The warp drive: hyper-fast travel within general relativity/ 
Class. Quant. Grav. 11 , L73-L77 (1994). 
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Appealing Characteristics 


Flat space-time inside the bubble 


Proper acceleration 
in the bubble is 
formally zero 


Images courtesy NASA 





Unappealing 
^ characteristic 

(square peg, round hole) 

MCC clocks synchronized with onboard 

clocks 


(divergence of phi = 0) 


(Coordinate time = proper time) 


Bubble Topology Optimization 

York Time magnitude decreases 




''bubble'' thickness decreases 


Energy density magnitude decreases 



' "bubble" thickness decreases 

10 

Surface plots of York Time & T°°, <v>=10c, 10 meter diameter volume, variable warp "bubble"" thickness 





As bubble thickness increases^ 
York Time intensity decreases 


Allowing the bubble to get 
thicker reduces the flat space- 
time real-estate in the center 


Changing topology greatly 
reduces the energy required 



But space-time is reaiiy stiff: c*/8tiG 


Can we further reduce the energy 
required by reducing the stiffness? 

Mavbe... but we need to engage 
higher dimensionai models to do so 
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Bubble Oscillation Optimization 








dt 
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♦ 
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dd) dU 

' 

dt dt 


Oscillate the bubble intensity 




Exotic Mass Warp Requirements, 10m diameter, 


Images courtesy NASA 










Warp Field Interferometer 






# 


He-Ne Laser 
~ 633 nm 


beam splitter 


7 ^ 


detector 



toroid capacitor ring 
normal to beam path 

/ resultant spherical field 



Warp Field Interferometer developed after 
putting metric into canonical form^: 





dt- 


Vsfir,) 


y;firS--i 


dx 


2 2 ^ 

~ dx + dy + dz 


Generate microscopic warp bubble that 
perturbs optical index by 1 part in 
10,000,000 

Induce relative phase shift between split 
beams that should be detectable. 


M(x, y)e 


2D Analytic 
Signal 
processing 


Numerical simulation Numerical simulation 


of fringe at detector of fringe at detector 
with device off with device on 

White, H., Discussion on space-time metric engineering," Gen. Rel. Grav. 35, 2025-2033 (2003). 




Numewc^J^m^ 
of contour pTot^f (j) 
for warp field 




Interferometer and Test-article Setup 






2 Polarizers 

(jointly used for 

intensity 

control) 


15 



PN 

NASA 


Fabry-Perot 

Interferometer 




Example: Michelson-Morley 
Interferometer image for 
Sodium source 





I 





Consists of two reflecting, 
highly parallel surfaces, called 
an Etalon 

The interference pattern is 
created within the Etalon 

Multiple reflections in the 
Etalon reinforce the areas 
where constructive and 
destructive interference 
occurs 

Allows for much higher- 
precision measurements of 
fringes (image averaging 
without software) 


Example: Fabry-Perot 
Interferometer image for 
Sodium source (note doublet) 





Frequency Spectrum (363,369) 



Frequency (Hz) 





Ml 



1§:01 


Primary/Auxiliary 


delator pressure 


Height control pressure 
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Open-air etalon Implementation 




Inactive 


Active 






f 


•\ 





Time of Flight Schematic 



200MHz test device 



Detectors 


Thorlabs High-Speed Avalanche Detector 



Model 

APD210 

Rise Time 

0.5 ns 

Supply Voltage 

+12 to +15 V 

Current Consumption 

200 mA 

Max. Incident Power 

10 mW 

Spectral Range 

400-1000 nm 

Frequency Range 

1-1600 mHz 

Maximum Gain 

2.5x10^ V/W 


HeNe 


EO Modulator 

Newport Broadband Amplitude Modulator 


Model 

Type 

Operating Frequency 


4102NF 

Broadband Amplitude Modulator 
DC-200 MHz 


Wavelength Range 
Material 
Maximum Vti 


500-900 nm 
MgOiLiNbOg 
195 V (5) 633 nm 


Maximum Input Power 2 W/mm^ (5) 532 nm 


Aperture Diameter 

2 mm 

RF Bandwidth 

200 MHz 

RF Connector 

SMA 

Input Impedance 

10 pF 

Maximum RF Power 

low 

Connector 

SMA 




Agilent Technologies Infiniium DS09254A 2.5 GHz Oscilloscope 

•2.5 GHz bandwidth across all 4 analog channels 
•20 GSa/s max. sample rate 

•Standard 20 Mpts memory per channel, upgradeable to 1 Gpts 
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LApiUlC LIIC Ul/yUL ucpci lUCI luy III lULUIC: lcdl 

devices 

- The idea of an optimized space warp needs vacuum 
energy, and large d^/dt 

— both of these conditions are predicted to be present in 
the q-thruster technology also being explored in the 


TMOlO 








Original Matthew Jeffries concept from mid 
1960% rendered by Mark Rademaker 


j -fT 



RademakerXdrex hies 


Matthew Jeffries is the artist that created the familiar Star Trek enterprise 
look 
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Act II: 

Experimental Thrust Measurement 

of RE Test Articles 
exploring q-thruster models 



Q-Thruster Background 

> A Q-thruster is a form of electric propulsion 

> Through the use of electric and magnetic fields, a Q-thruster pushes quantum 
particles (electrons/positrons) in one direction, while the Q-thruster recoils to 
conserve momentum 

> This principle is similar to how a submarine uses its propeller to push water 
in one direction, while the submarine recoils to conserve momentum 

> Based on test and theoretical model development, expected thrust to 
power for initial flight applications is 0.4N/kW (~7x Hall) 

> 0.4 N/kW enables power-constrained HEO SEP missions to close without 
needing chemical kick stages and very long transit times. 

> 0.4N/kW coupled with persistent power (e.g. NEP) enables rapid transit 
missions throughout the solar system. 


'Tf vacuum fluctuations can be harnessed for propulsion by anyone besides 
science-fiction writers, the purely engineering problems of interstellar flight 
would be solved." - A.C. Clarke 
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Cannae RF Test Article 

> Test article is a pillbox/beam pipe design 
fashioned after an RF resonant cavity design 
used in high energy particle accelerators. 

> Each test article is ~11" in diameter and 4-5" 
between ends of the beam pipes, not counting 
beam pipe extensions or antenna mounts. 

> Concept is a modified particle accelerator 
design that incorporates engraved radial slots 
along the outside edge of the resonant cavity 
interior, but on only one side of the pillbox 
(equatorially-asymmetric) . 

> Cannae theorizes that asymmetric engraved 
slots would result in a force imbalance 
(thrust). 

> As a result, a second (control) test article was 
fabricated without the internal slotting 
dubbed the null test article. 
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5.0V Xmit 
On Off 


0-to-28Vdc 
Freq. Control 


O-to-lTVdc 
Atten. Control 


28Vdcat O-to-11.0 
Amps Max 


VOM#l VOM#2 


(Ten?>) 


(2. 22AWG) 


(2. 22AWG) 


ZX95 1015-S+ 
750-to-l,010 MHz 
0-to-28Vdc VCO 


ZX73-2500M-S^ 
» lO-to-2,500 MHz 
W-Atten. 


ZHL-100W-13+ 
800-to-l,000 MHz 
lOOW (50dBm) Amp 


ZX47-40LN-S+ 
10-8,000 MHz 
Power Detector # 1 


Cannae Fwd. / Rev. 
Directional 
Coupler DRl 


Spectrum 

Anatyzer 


ZX47-40LN-S+ 
10-8,000 MHz 
Power Detector #2 


■Whip 

Antennas 


^28V-to-5.0Vdc 


CivST ff. LM Contacts & #22 
AWG wiring implies 7.0 Amps 
max per chassis circuit. 


VOM#4 i (No\ 


VOM#3 


Testing 

Overview 


thrust 


thrust 


RF signal fed into one of test 
article beam pipes via an 
approximately 20-gauge copper 
wire drive antenna. 

On opposite beam pipe, an 
approximately 20-gauge copper 
wire serves as a sense antenna 
for manual tuning during 
testing. 

The longer beam pipe is the RF drive 
antenna that in practice ends up being 
a ^4 wave resonance system in its own 
right and has a dielectric PTFE slug in 
the throat in both the slotted and null 
test article. 

It is this characteristic that became an 
item of further consideration after 
completion of the test campaign. 


Distance and Chart Data 


Thrust Pulse = ~33 uN 
i (Up = to Left) 


Thrust 


jPulse ON Time 


200V Cal Pulses =~26 uN 


(Down = to Right) 


196 Seconds (5.1 Points/sec at 2048 Samples/Point) 
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L - 71.9 



Test article does not get 
warm during operation 
(high Q resonance system) 
Amplifier kept below 100 F 


Cannae Test Results Summary 

> The resistive RF Load evaluation indicated no significant 
systemic cause for torsion pendulum displacement. 

> Based upon this observation, both test articles (slotted and 
unslotted) produced significant thrust in both orientations 
(forward and reverse). 

> Test schedule constraints prevented multiple data points 
to be gathered in the reverse orientation, and the single 
data point for each test article is insufficient to allow 
comparative conclusions (between slotted and unslotted) 
to be drawn. 

> However, for the forward thrust orientation, the difference 
in mean thrust between the slotted and unslotted was less 
than two percent. 

> Thrust production did not appear to be dependent on the 
slotting. 



Configuration lA 



Configuration 

Test Article 

Thrust 

Direction 

Thrust Range 
(fxN) 

Mean 

Thrust 

(uN) 

Number of 
Test Runs 

lA 

Slotted 

Forward 

31.7-45.3 

40.0 

5 

IB 

Slotted 

Reverse 

48.5 

48.5 

1 

2A 

Unslotted 

Forward 

35.3 - 50.1 

40.7 

4 

2B 

Unslotted 

Reverse 

22.5 

22.5 

1 

RF Load 

50Q Load 

N/A 

0.00 

0.00 

2 



Configuration 2B 


RF Dummy Load 
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COMSOL Analysis of Cannae Test Articles 


> 


Computer modeling of the electric field 
within the pillbox and beam pipe 
illustrates weakness of electric field in 
vicinity of cavity slots and relative strength 
of the electric field within the beam pipe, 
especially in the drive antenna coaxial 
cable and the region around the cable 
within the PFTE dielectric slug as seen in 
the following figure. 



> Consideration of the dynamic fields in the k4 wave resonance tube shows that there 
is always a net Poynting vector meaning that the RF launcher tube assembly with 
dielectric cylinder common to both the slotted and smooth test articles was 
potentially a Q-thruster where the pillbox is simply a matching network. 


> The predicted thrust using the Q-thruster analytic model is 34 micronewtons with 
input power of 25 Watts and Quality factor of 8000. 
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Tapered RF Q-thruster 

> Just like there are many ways to build a race car motor, there are 
many ways to build a thruster, or in this case, a Q-thruster. 

> The figures depict one of the early COMSOL models representing 
an early possible construction of a tapered RF unit alongside the 
actual construction that was finally implemented as informed by 
the COMSOL findings. 

> The tapered RF geometry was explored as it pertains to some recent 
findings published by the Northwest Poly technical University. 



Test article 
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4" dielectric 
resonator 


TE012 mode 
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I 

Browse i 
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577,0001 

576*000 

575.000 






Points Saved 0 


File Name 
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COMSOL RF Analysis 
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TM212, 1,937.1 15MHz, ~35W, Ql=6,726, 


Microns Distance and Chart Data 
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TM212, 1,937.1 15MHz, ~50W, Ql=6,726, 
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Forward Thrust Campaign 





FORWARD WORK 

> JSC: Complete the vacuum test campaign 


> Perform an IV&V test campaign at the Glenn Research Center using 
their low thrust torsion pendulum followed by a possible repeat 
campaign at the Jet Propulsion Laboratory using their low thrust 
torsion pendulum. 


> The Johns Hopkins University Applied Physics Laboratory has also 
expressed an interest in performing a Cavendish Balance style test 
with the IV&V shipset. 
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GRC Torsion Pendulum 




Act III: 

Value Proposition 



■ Dawn propulsion & 
power system 

■ Three NSTAR xenon 
ion thrusters 

■ 0.04 N/kWe thrust- 
to-power ratio 

■ ~90 mN/thruster 

■ 2.3 kWe solar power 
atl AU 

■ 425 kg of xenon 
propellant 

Dawn with 2.3kWe Ion Thrusters Dawn with 2.3kWe Q-thmsters 



Activity 

Dates 

Segment 

(months 

) 

Elapsed 

(months) 

Dates 

Segment 

(month 

Elapsed 

(months) 

Earth to 
Vesta 

Sep 2007- 
Jul 2011 

46 

46 

Sep 2008 - 
Apr 2009 

7 

(15% of ion) 

7 

Vesta orbit 

Jul 2011 - 
Jul 2012 

12 

(fixed) 

58 

Apr 2009 - 
Apr 2010 

12 

(fixed) 

19 

Vesta to 
Ceres 

Jul 2012- 
Feb 2015 

31 

89 

Apr 2010 - 
Jul 2011 

15 

(48% of ion) 

34 

Ceres orbit 

Feb-Jul 2015 

5 

(fixed) 

94 

(7.8 yrs) 

Jul-Dec 2011 

5 

(fixed) 

39 

(3.3 yrs) (41%) 
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300 kW SEP Mars 


• 70t stack departs from 
DRO 

• 300kW SEP 

• 0.4N/kW Q-thrusters 

• 50-day stay in Deimos 
orbit around Mars 

• Total mission duration of 
788 days and 2AU 
maximum distance from 
the sun. 
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2MW NEP Mars 


J^ARS 

66 days 

ARTH 


UN 


• 90t spacecraft 

• 2MW power 

; 0.4N/kW (SOON) 

• 246 day mission with 
70 day stay at Mars 



50t cargo 


• 20t 2MW nuclear reactor (10 kg/kW) 


• 20t Q-thruster bank (lOkg/kW) 
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Crewed Titan/Enceladus Mission 

• Same vehicle characteristics as Jupiter mission assumed 

• Trip time to/from Saturn ~9 months each way, 6 month stay in Europa 
vicinity, 6 month stay in Enceladus vicinity. 32-month total mission. 

• Shorter than "current" conjunction-class Mars missions 



Assuming: 

J=0A N/kWe 
P=2000 kWe 

«power=10 kg/kWe 

ap,op=10 kg/kWe 
nfipi = 50t 

T/m =0.91 m-g's 
Vehicle mass = 90 1 




Act IV: 

Dynamics of the Vacuum and Casimir 
Analogs to the Hydrogen Atom 



Principles of Q-thruster Operation 


> 


> 


> 


Principle 1 : Local mass concentrations, say in the 
form of a conventional capacitor with a ceramic 
dielectric, affect vacuum fluctuation density 
according to equation 1 


Principle 2 : Just as relativistic acceleration (LFnruh 
radiation) can change the apparent relative density 
of the vacuum, so too can higher order derivatives 
according to equation 2. 


Principle 3 : The tools of MagnetoHydroDynamics 
(MHD) can be used to model this modified 
vacuum fluctuation density analogous to how 
conventional forms of electric propulsion model 
propellant behavior. 
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TABLE I: This table shows the derived '^density’’ of a 


given energy state n, with Z = 1 


71^ 

7'adius(m) 

E{eV) 


E{J) 


p {kg/ 

m^; 

) 

1 

5. 

29 

X 

10“ 

TT 

13.60 

2, 

.176 

X 

10" 

-18 

3.905 

X 

10' 

-5 

2 

2. 

11 

X 

10“ 

■10 

3.40 

5, 

.440 

X 

10" 

-19 

1.525 

X 

10' 

-7 

3 

4. 

76 

X 

10“ 

10 

1.51 

2, 

.418 

X 

10“ 

-19 

5.952 

X 

10' 

-9 

4 

8, 

46 

X 

10“ 

■10 

0.85 

1, 

.360 

X 

10“ 

-19 

5.959 

X 

10“ 

10 

5 

1 

.32 

X 

10' 

-9 

0.54 

8, 

.704 

X 

10“ 

-20 

9.997 

X 

10“ 

11 

6 

1 

.90 

X 

10' 

-9 

0.38 

6, 

.044 

X 

10“ 

-20 

2.325 

X 

10“ 

11 

7 

2 

.59 

X 

10' 

-9 

0.28 

4, 

.441 

X 

10“ 

-20 

6.774 

X 

10“ 

■12 


^ The primary quantum number n is only varied from 1 to 7 here. 


47T€on^h^ 

e^irie 

= 5.29xl0“^^7ne^ers, n = 1,2,3,.., 


En — — 


1 


m 

( VI 

1 


V 47760/ 



= ^1.3.6eV, 1,2,3, 




^ Zjji 



5.29x10 ^bneters, 77 = 1,2,3,.., 



^^13.6el/, n = 1,2, 3, ... 




n 


9 4 

^ 2^777 


.3 

Z,n 
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TABLE II: This table compares the derived “density” 
of a given energy state n, with Z = 1 to the Casimir 
density for a cavity with a separation distance of 2r„. 


n 

radius{7n) 

p {kg/m^ 

) Casivnir {kgj 


Ratio^ 

1 

5.29 

X 

10“ 


3.91 

X 

10' 


1.16 

X 

10' 

-4 

2.96 

2 

2.11 

X 

10“ 

10 

1.53 

X 

10' 

-7 

4.51 

X 

10“ 

-7 

2.96 

3 

4.76 

X 

10“ 

10 

5.95 

X 

10' 

-9 

1.76 

X 

10“ 

-8 

2.96 

4 

8.46 

X 

10“ 

10 

5.96 

X 

10“ 

-10 

1.76 

X 

10- 

-9 

2.96 

5 

1.32 

X 

10' 

-9 

1.00 

X 

10“ 

-10 

2.96 

X 

10“ 

■10 

2.96 

6 

1.90 

X 

10' 

-9 

2.33 

X 

10“ 

-11 

6.88 

X 

10“ 

■11 

2.96 

7 

2.59 

X 

10' 

-9 

6.77 

X 

10“ 

-12 

2.00 

X 

10“ 

■11 

2.96 


^ This is the Casimir force per unit area multiplied by 1/c^. 

^ This ratio is the Casimir column value divided by the p column 
\'alue 


Friedmann equation: 


d 

a 


47tG 


3 


{pc^ + 3P) 


P < -p(?j% 
w — 1 / 3 




3 

Z,n 


1 hcw^ 


ip) 


3c2 240d'‘ 



Vacuum density perturbation (kg/m'^3) 


Vacuum Perturbation 



l.OOE-12 l.OOE-11 l.OOE-10 l.OOE-09 l.OOE-08 

Radius (m) 


P [kg 1 771^) 



radius {meters) 


-10 


-9 



-0.5 


-1.5 


-2.5 


Mesh sizes: 
n=l: 1 pm 
n=2: 2.5 pm 
n=3: 5 pm 
n=4: 10 pm 
n=5: 25 pm 
n=6: 50 pm 
n=7: 50 pm 


-2-101234 


FIG. 6: COMSOL 2D axisymmetric model: element 
sizes for n = 1 to = 7 is Iprn, 2.5pm, 5pm, 10pm, 
25p77i, 50p7?7, and 50p/7i respectively. In this figure, the 
mesh size is too dense to be discernable for 77 = 5 and 
lower. FIG. 4 shows the mesh for the 77 = 1 and n = 2 


regions 



FIG. 4: Close-up of COMSOL 2D axisymmetric model. 
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(b) Microview 


FIG. 5: COMSOL analysis results for 77 = 1 
eigenfrequency: panel 5a shows the model out to the 
= 6 orbital, and panel 5b shows a close-up view' of the 
77 = 1 solution 





TABLE r\^: This table shows the predicted “orbital” 
frequency and the COMSOL analysis eigenfrequencies 
for the n = 1 to n = 7 states for hydrogen. 


n 

Orhital f req^ 

COMSOL freq^ 

% err or 

1 

6.. 58 

X 

10^^ 

6.25 

X 

10^^ 

-4.98 

2 

8.23 

X 

10^^ 

8.23 

X 

10^^ 

0.05 

3 

2.44 

X 

lO''* 

2.38 

X 

10'^ 

-2.48 

4 

1.03 

X 

lO'^* 

1.01 

X 

10^^ 

-1.59 

5 

S.26 

X 

10^^ 

4.98 

X 

10^^ 

-.5,36 

6 

3.05 

X 

10^3 

3.48 

X 

10^3 

14.28 

7 

1.92 

X 

10^^ 

2.13 

X 

10^^ 

11.16 


^ Orbital frequency is in Hz. 

^ COMSOL frequency is in Hz. 
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FIG. 7: COAISOL analysis results of the acoustic 
‘Tiaturar’ vaciiuin model The orbital shells (dark lines) 
can be counted, but the n = 1 radius is quite small as 
seen from the top left thumbnail that depicts the 
COMSOL eigeufrequeiicy solution for that orbital. 
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FIG. 9: 2D Axisymmetric model results that capture 
axisymmetric acoustic solutions like 2p, 3d, and 4/ 

orbitals (m = 0). 



2 m dx^ 


-|- V Ip 



FIG. 3: Plot for the Z = 1, 2p orbital from Orbital 
Viewer Software. 





